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, Abstract 
CN ■ 

We use simple iterated one-loop graphs in massless Yukawa theory and QED to pose 
psj , the following question: what are the symmetries of the residues of a graph under a per- 

^ ' mutation of places to insert subdivergences. The investigation confirms partial invariance 

of the residue under such permutations: the highest weight transcendental is invariant 
under such a permutation. For QED this result is gauge invariant, i.e. the permutation 
invariance holds for any gauge. Computations are done making use of the Hopf algebra 
structure of graphs and employing GiNaC to automate the calculations. 



1 Introduction 



Q^l This paper serves three purposes: i) it employs GiNaC in Feynman diagram calculations 



and provides algorithms which automate the renormahzation process, very much in the spirit 
of Q; ii) it investigates symmetries of short distance singularities under permutations of places 
where to insert subdivergences in a graph; iii) it once more confirms the presence or absence 
^ ■ of transcendental coefficients of short-distance singularity in accordance with the topology of 

^ . a graph. 

Our laboratory of investigation are simple one-loop Feynman graphs in massless Yukawa 
theory or QED, inserted into each other in tree-like hierarchies. Thus, the combinatorics of 
renormahzation boils down to the Hopf algebra of decorated rooted trees with only a small 
number of decorations and the analytical challenge posed by those decorations reduces to 
expansions of F- functions near unit argument. The question we can ask is for the distribution 
of the Riemann C-function over the various poles in the Laurent series of graphs of that form. 

In contrast to Q and subsequent papers Q, where the renormahzation problem was 
automated in a similar context optimized for speed and efficiency, we have developed here 
algorithms which allow for non-trivial spin-structures and an easy generalization to arbitrary 
decorations: the primitive decorations can be inserted as arguments so that the algorithm 
can handle arbitrary primitive graphs when their analytic structure becomes known. 

We work in the context of dimensional regularization, so that any Feynman graph becomes 
a Laurent series in e = (4 — D)/2, the deviation of the dimension from its integer value four, 
and the pole terms reflect the short distance singularities in the theory. The flrst order pole is 
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denoted as the residue of a graph. Its significance hes in the fact that higher pole terms can 
be reduced to polynomials in residues and in the fact that the residue of a primitive graph 
is an invariant under diffeomorphisms of external parameters of the graph (diffeomorphisms 
of external momenta and masses) as well as an invariant under variations of renormalization 
schemes. Such a residue is a motivic number then in some modern mathematic parlance. The 
question as to which class of such numbers is sufficient to describe the residues of a given 
quantum field theory is open and fascinating Q. 

We only study two much more basic questions, motivated by previous and ongoing inves- 
tigations into the analytic structure of pole terms and residues in particular. 

The first is the independence of the appearance of transcendentals under variations of the 
quantum field theory which realizes a graph with a given topology. To specify the topology 
of a one-particle irreducible Feynman graph F, let us consider the adjacency matrix M(F) of 
F. If F has n vertices, this matrix is a n x n matrix. We take for its non-zero entries pairs 
(propagator type, powercounting weight), ie. each non-zero entry M{T)ij specifies that vertex 
i is connected to vertex j by a propagator of some type, which has a certain powercounting 
weight .0 In the cases studied here, the possible entries are 

(fermion,!), (photon, 2), (scalar boson, 2). 

The listing of the powercounting weight is redundant, as it is determined by the type of the 
propagator. We list it just for easy reference. 

The graphs of QED and Yukawa theory which we will compare always have adjacency 
matrices which agree in all their zero entries, 

M(FQED)ii = 0^ M(FYuk)ii = 

and agree for each non-zero entry in the powercounting degree of the corresponding edge. 
It is only the nature of the edges which changes from the spin one vector boson propagator 
-the QED photon- to the spin zero scalar boson propagator in Yukawa theory. Note that the 
structure of short distance singularity then remains fixed in the transition of one theory to 
the other. We then expect and confirm that rational numbers can vary in the transition from 
one theory to the other, while the transcendentals we see remain invariant and specific for the 
chosen topology. Here, a topology is uniquely described by considering in such an adjacency 
matrix all non-zero entries as equal. So it just gives information about how vertices are 
connected, but forgets about the nature of the propagators establishing that connection. To 
find non-rational numbers, we have to go up to four loops at least in the simple class of 
Feynman graphs which we consider. There, it is the swiss cheese topology of Fig. la in which 
we expect to see a residue ~ C(3), while the ladder topology of Fig. lb are known to have only 
rational residues [^. 

The second question is of different nature: for a log-divergent vertex graph which is prim- 
itive under the coproduct, it is evident that its residue has the above mentioned invariances. 
The Hopf algebra structure immediately allows to prove that in a vertex graph which is not 
primitive under the coproduct such that the graph contains divergent subgraphs, the coeffi- 
cient of the highest order pole still has these symmetries and is given by an easy calculable 
product of residues with a combinatorial factor, determined by the scattering type formula 
of Q, which incorporates the 't Hooft relations between higher pole terms in graphs. 

^We can notate the type of vertex in the diagonal entries of this symmetric matrix (no vertex is connected 
to itself by a propagator). 
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What we call here the 't Hooft relations in accordance with |Q is the simple fact that 
higher order poles in a graph T with UV-divergent subgraphs 7 can be calculated from prod- 
ucts of lower order poles of these subgraphs 7 and their complements r/7. This is well-known 
and a necessary requirement to make the renormalization group work: the Z-factor for a given 
physical quantity is an invertible (it starts with 1) formal series over counterterms of graphs 
such that its logarithmic derivative with respect to the variation of the renormalization scale 
u, d\og{Z) / du"^ , is finite. This finiteness establishes relations between pole terms which, for 
the case of dimensional regularization, were, it seems, first explored by 't Hooft. These rela- 
tions are a direct consequence of the mathematical structure of the Hopf algebra underlying 
renormalization, and its one-parameter group of automorphisms 

Let us now consider the residue of such a graph which does have divergent subgraphs. 
Typically, this residue will be a number which can be decomposed in terms of transcendental 
weight: it will contain contributions ranging from rational numbers to monomials in (,{]) 
of up to transcendental weight ? — 1, where / is the bidegree of the graph, calculated from 
its coproduct Q^. We cannot expect the whole residue to be invariant under the above 
symmetries, as rational numbers can and will vary freely. But here we report a remarkable 
partial symmetry observed in our rather restricted lab of iterated one- loop graphs: the highest 
weight transcendental in the residue is invariant under permutations of external momentum 
as described below. We confirm this by empirical calculation to high loop orders. We finally 
prove the result in the context of the simple iterated one-loop graphs considered here. One has 
an almost elementary proof in this context and we will discuss the difficulties which arise in 
the general case. The nature of this result fits nicely with a structural investigation of Dyson- 
Schwinger equations to be delivered elsewhere [^. Note again that changing the momentum 
flow but maintaining the topology corresponds to alterations of the type of non-zero entries in 
a suitable adjacency matrix, as Fig.§ clearly exhibits. Again, the degree of powercounting and 
hence the structure of short distance singularities, as well as the topology, remain unchanged. 

To understand the type of symmetry we want to investigate, let us consider the one-loop 
vertex function. If T{pi,p2,pz) is a vertex correction with momentum pi for the external 
boson, p2 for the incoming fermion and Pz = Pi + P2 for the outgoing fermion, then we will 
compare := T{0,p,p) with T^ := r{p,p,0) at the one-loop level. It is a single permutation 
of the flow of an external momentum for a vertex function at zero momentum transfer, 
which distinguishes the two one-loop functions. Locality of counterterms ensures that the 
counterterm for a vertex correction graph is invariant under this permutation for any graph. 
So nothing exciting can be learned from just changing the momentum flow in a given vertex 
graph. Now let us start with such one-loop vertex corrections and let us insert more one- 
loop vertex corrections always at the vertex of zero momentum transfer (zmt), which indeed 
changes under the above permutation. This gives the topological equivalent Feynman graphs 
F^,F^, with permuted types of propagators, of Fig.^. They indeed have adjacency matrices 
such that again 

M{r% = ^ M{r% = 0. 

Note that it is the requirement to maintain the same topology which forces us to consider 
quite different looking graphs. For them, the usual field theoretic requirements only demand 
that the coefficients of the highest order poles are identical. All non-leading pole terms in e 
do not have to and indeed do not agree. 

■^For all graphs considered here, it simply agrees with the number of loops in Yukawa theory. For the QED 
vertex, a shift by one unit can appear as explained later in the text. 
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(a) "Swiss cheese" topology (b) "ladder" topology 

Figure 1: "Swiss cheese" and "ladder" master topology. 

But a closer investigation will establish a remaining partial symmetry between those 
graphs, such that the highest weight transcendental in the residue is invariant. Such sym- 
metries are of great interest. Feynman graphs with sub divergences can be built by using the 
pre-Lie structure of graphs, which results from underlying insertion operads Q. The question 
to what extent the permutation group acts trivially under such insertions is a natural question 
in operad theory which needs to be answered to understand these operads. It directly leads 
to the questions studied here in a simplified context. We regard the results reported here as 
a first step in a detailed analysis of actions of the permutation group in this context. 

With regard to the first question, we remind the reader that for a massless Yukawa theory 
represented as a Laurent-series in e, it was already shown in |^] that graphs with the swiss 
cheese topology of figure [l|.a possess a Cl'^) in their counterterm, whereas graphs with a ladder 
topology like in figure [l|.b just evaluate to rational coefficients. Similar results were obtained 
in QED and other theories for vertex functions of type T{0,p,p) and self-energies, in [^, ^, 

Our first aim is now to find out if this remains true in massless QED for both types 
of momentum fiow and if there are any symmetries in the coefficients of the ^-functions 
in both theories, under the permutation between the two types of vertex functions and 
r^. Therefore we will rebuild the scheme given in |Q for the massless Yukawa theory with an 
extension to vertex corrections that carry a different fiow of momentum in the sense described 
above, and similarly for QED using the matrix calculus of |]lO|. The resulting algorithm is 
implemented using the GiNaC library and will be described in section ^. With the help of 
this program, we will calculate the antipode of graphs that represent the different topologies 
given above and compare them. 



2 Calculations 

Consider a one-loop contribution to the fermionic propagator: 



where the subscript "[0,0]" will soon serve to give the number of one-loop sub divergences to 
be inserted inside the diagram at the fermionic or bosonic line. All integrals considered in 
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this paper can be reduced to integrals Fa^ 

I{q;a,b) = -^[d^'k 



_l_[,2](2-(a+.)-e)^^^^ (1) 



with 



_ r(2-Q-£)r(2-6-£)r(Q + 6-2 + £) 

ra,b~rb,a - r (a)r (6)r (4 - a - 6 - 2e) ' ^ ' 

which we typically need for a = ni + n2£,b = mi + m2S, for integers ni, ?i2, mi, m2. This 
reduces the identification of non-rationals (transcendentals, we dare say in the following) to 
an expansion of the F-function near unit argument, as promised. Note that 

Fa,-n = i^-n,6 = , nGN°. (3) 

Accordingly, in our conventions 

^o,o]{q') = [qT'Ifi^i i =■■ [g']-^So,o ^ (4) 
2.1 Yukawa theory 

The study of iterated one-loop integrals reduces to the study of the following elementary 
functions, which we will call characterizing functions: The one-loop fermion self-energy with 
insertions at its two internal propagators demands knowledge of 

Si J := -[Fj£^i+j£ -I- Fi+ie^i+je — Fi+jeje]. (5) 

The one-loop boson self-energy with insertions at its two internal edges is given by Hij and 
the two one-loop vertex functions at zero momentum transfer with insertions at either the 
zero momentum vertex, or internal lines, need knowledge of functions Tij for the vertex 
corrections. 

nij/[tr(l)] := -l[Fi,,i+,-,-Fi+i,,i+,-, + Fi+i,,,,] (6) 

■= -^l+ie,l+i£ (7) 
2 1 

'■= 2^F"2+ie,je - F2+ie,l+je + Fi+ie,l+je]- (8) 

Here, we have divided the boson self-energy by the trace of the unit matrix tr(l) (trace over 
spinorial indices) for easier comparison of insertions of subgraphs into bosonic and fermionic 
lines. 



■^These F-functions will give rise to the ("-functions we are looking for. An easy way to see the connection 
between the F- and the (^-function is the formula: 

F(l + e) = exp(-7e) exp (^f;^(-£)"j , le| < 1 

Internally, GiNaC follows a different approach. It computes the derivatives r'(a;) = r{x)ip(x) and -^j-^tpix) = 
ipnix) in terms of polygamma functions tpmix) and "knows" how to evaluate polygamma functions at integer 
arguments, e.g. — —7 and V'n(l) ~ {—)"'^^n\(^{n + 1). 
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(a) Vertex correction of type one. 



q 

(b) Vertex correction of type two. 



Figure 2: The two types of vertex corrections. 



The functions and F'^ represent the two types of vertex corrections shown in figure |2|. 
The indices i and j give the number of subdivergences at different Hues: 
Sjj-: i = number of subdivergences at the fermion hne, 

j = number of subdivergences at the boson hne. 
Hi J- i = number of subdivergences at the lower fermion hne, 

j = number of subdivergences at the upper fermion Hne. 
F^^: i = number of fermion self-energies and vertex corrections 

plugged into the zmt vertex and the internal edges connected to it, 

j = number of subdivergences at the boson line not connected to this vertex. 
F?^: i = number of fermion or boson self-energies and vertex corrections 

plugged into the zmt vertex and the internal edges connected to it, 

j = number of fermion self-energies at the fermion line not connected to 

the vertex of zmt. 

To denote complete graphs, we introduce additional functions ^[ij], H[i,j] and ^[ij^k] which 
notate the different kinds of subdivergences (cf. ^[ofi] in @). The indices count for ^[ij] and 
Iljj j] as the indices in their corresponding characterizing functions given above. For vertex 
corrections of type one and two, the notation Fjj ^^^j denotes the three types of insertions in 
the following way: i insertions of self-energies at the edges connected to the zero momentum 
vertex, j insertions of self-energies at the edge not connected to that vertex, and k insertions 
of vertex corrections at the zmt vertex. Table || gives examples. 



2.2 QED 

The main difference stem from the presence of a gauge parameter ^ in the photon propagator 
and from the — ie7^ vertex which make the calculations more difficult. In fact, as we will see. 



we will have to deal with matrices for the vertex corrections [10|. Nevertheless, the structure 
for the translation of a graph to an analytical result will be similar to what we had in Yukawa 
theory. Hence we will mainly give the results and just list the relevant changes. 



2.2.1 Vacuum polarization 

We will only consider one-loop photon self-energies, vacuum polarizations, as in our restricted 
class of Feynman graphs we can not construct a gauge invariant set of vacuum polarizations 
at higher loop orders. 



6 



-o- 




^[0,01(9') 

S[1,0]('?') 

n[o,o](g') 
n[2,i](9') 
^[0,0,0] (9^) 




^f0.0.2l(9^) 








[0,0,2] » 
^[2,0,01(5^) 

to,i,o]( 



rfo,i,oi('?') 



^fo,o,o](9^) 

^fo,0,2](5^) 

rfi,i,i]('/^) 



-ig 
-ig 
-ig 
-ig 
-ig 

-ig 

-ig 

-ig 

-ig 

-ig 

-ig 



8s2,ino,o(So,o)2[g2]-^^^ 
8n2,iSi,o(So,o)'[g']-V 

7 pi pi pi [„2i-3e 
2,0-'- 1,0-*- 0,0W J 



^2,0^1,0^^0,0 [9^] 



21 -Se 



5 -nl 



21 -2£ 



rJ.ino,ok1 



7p2 p2 p2 r 
^ 2.0-*- 1,0-*- 0,0 W 



21 -3e 



'rinrL(So,o)^['z 



21 -4e 



Table 1: Some examples how self-energy and vacuum polarization graphs are built up in Yukawa 
theory. 

2.2.2 Fermion self-energy 

The fermionic propagator in QED with insertions at the fermionic and photonic lines demands 
E: 

^ 1 



J{2-D) -{3-D) Fi+j,^i, - (3 - 

+ F2+j£-l+ie — 2 F2+je,i£ + -^2+je,l+ie ] (1 ~ '^0,j) 

+ [(2-I))Si,o + eS-o]<5o,i, 



with 



1 r 

- [i^l+ie,l + Fie,l] 



2 [-^-l+ie,2 — 2Fje,2 — Fie,! + Fl+ie,2 — ^l+ie,l\ 



F 



(9) 

(10) 
(11) 



The appearance of the Kronecker 5o,j is obvious from the fact that the presence of one-loop 
vacuum polarizations in the internal photon line forces transversality of that propagation. 
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Figure 3: One- loop contribution to the vertex correction. 



2.2.3 Vertex corrections 

The most important difference to tlie Yukawa theory occurs in this part of the calculation. 
Consider the one-loop vertex correction in QED for vertex graphs as shown in Fig. |3|. 

The external structure of the vertex consists of two form factors: One for 7^ and one for 
as T^{0,q,q) = Fi{q^)^^ + -^2(9^)^- In the previous calculations, every subdivergence 
just caused an change in the exponent of various momenta in denominators of the integrals, 
resulting in non-integral exponents in equation (||). Now the result has two form factors Fi, F2 
that emerge through the evaluation of a subgraph in a "bigger" graph. We will administer 



the two form factors in a matrix notation in accordance with |1C]. We will need two-by-two 



matrices. The four entries determine four functions A^*'^^ Here, a, 6 € {1,2} and count 
the number of internal insertions as before. The case 6 = 1 corresponds to an internal vertex 
7^, the case 6 = 2 corresponds to an internal vertex q^xi/q^ (where q, say, is the momentum 
flowing through this zero- momentum transfer vertex), while the index a enumerates the two 
possible form factors in the result. The result for the one-loop graph of Fig.^ is then 



^(0,0) , ^(0,0)_ 
^1,1 7At ^ ^2,1 ^2 



For the two-loop graph shown in figure we find 



[0,0,1] 



. (0,0) , (1,0) ^(0,0) . (1,0)\ / . (0,0) , (1,0) A (0,0) . 

^11 ^11 ^ ^21 ^12 y 7m ^21 + ^21 ^ 



(0,0)^(1,0) 
22 



-2]-2^ (12) 




i q-k 



Figure 4: Two-loop contribution to the vertex correction. 



Hi) 

The multiplication of the A^^ can be reformulated as a matrix multiplication. We define 
the general matrix (the upper index refers to the two cases of vertex corrections which we 
consider, we omit that index for simplicity in the matrix entries): 

"^'^ - I, A^^) A^^ ) 

To each part of a graph which has the form of a vertex correction we assign a matrix. In 
this case of vertex corrections of type one, as in ^.a, the index i is the total number of 
sub divergences at the fermionic line (with no difference if it is of the form T, or T) and j 
the number of subdivergences at the photon line, as in Yukawa theory. One multiplies the 
matrices for the different divergent parts of a graph, starting from the innermost vertex 
correction which contains the vertex of zmt, and the corresponding matrix has zero entries in 
the second column obviously. 

Let us make this clear with the help of our example of the two-loop vertex correction 
again. For the graph of figure ^ this means that we begin with the inner vertex correction 
marked with a box: 




and then multiply this with the matrix for the outer vertex correction, which has one vertex 
correction as a sub divergence: 




/aS^) A^^M/Af/) 0\_/ Ail'°)AS°) + AS°)AS'°^ 
^ 1^ A^f^ a£°) ) [ aT o)-[ A^^f^AT + A^'^Aff^ 

Muhiplying (||) with (1, 0)^[g2]-2e, we get the result as the column vector 



(13) 



/(Asi-°)Ar+A&°)A(r)[,r^M 

[{A^^^^AT+A^^^ATm-'^ )' 

where in general the "upper" entry of this vector is the form factor Fi(g^) belonging to 7^ 
and the "lower" one the form factor F2{q^) belonging to 

Subdivergences that are not vertex corrections given by Sj^j and Xlj j are multiplied with a 
unit matrix and inserted in the string of matrices in front of the matrix of the vertex correction 
which they are part of. They increase subscripts i,j accordingly in that matrix. 

As an example we get for the following graph: 
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A (1,0) A (1,0) 

^11 ^12 

A (1,0) A (1,0) 

^21 ^22 

= M2_;^no,oM^^,o^o,o- 

As already mentioned, the indices i and j count for QED vertex corrections of type one, 
^[ijfc]' similar to the ones defined for Yukawa vertex corrections of this type. Similarly, the 
indices for the type two vertex corrections r| j ^ in QED count as in Yukawa type two. And 
the matrix multiplication described above stays the same apart from different entries for the 
matrices, which can be readily calculated in terms of functions -Fa, 6- We do not list the two 
sets of four functions A^*'J^ in terms of Fc,d functions explicitly, but the interested reader 



can find them from our publicly available programs H. Table I shows some examples how 
graphs are built in QED. 



3 Renormalization 

Renormalization employs a simple principle of multiplicative subtraction, making use of the 
underlying Hopf algebra structure of Feynman graphs [p!^ , |^ : the coproduct 

A(r) = r0i + ior + ^7or/7 

and antipode 

5(r) = -r-^7r/7 

7cr 

are the structure maps which allow the construction of counterterms and renormalized quan- 
tities. One employs Feynman rules <j) : H ^ V as an element in the group of characters 
of the Hopf algebra H, with target space V (a suitable ring or algebra) and makes the 
target space into a Baxter algebra {V, R) by choosing a renormalization map R such that 
R{ab) + R{a)R{b) = R{aR{b)) + R{R{a)b). One then has the counterterm 

5^(r) = -R[m + sii^mh)] 

7cr 

and a further recursion delivers the renormalized result 

si * m = 4(r) + m + ^ sJ(7)0(r/7). 

7cr 



10 



^[0,0] 19 ) - 






^[1,0] w J — 


(-ie)^Si,oSo,o[g^ 




■^[2,0] W ) — 


{-ieft2,oto,oto,o[q^r^'^ 


■^[2,1] W ) — 


(-ie)8S2,ino,o(S 




■'-'-[0,0] w ) — 


(-ie)2[g2]-eno,o 




[0,0,0] w J - 


(-ie)3Mio[9']-^ 




[0,0,2] W j - 


i-ief MloMl,Ml,[q'']-^^ 


[2,0,0] W J - 


{-iey Ml^^ti^o^ 


o,o[gr'^ 


r[o,i,o]('?') = 




^fo,o,o](9^) = 


{-igf Ml,[q^]-^ 




^fo,0,2](9^) = 


i-igY Ml,Ml,Ml,[q^]-^^ 


rfi.i,i]('^^) = 


i-ieyMl^Mlitofinq']-"' 



Table 2: Some examples how self-energy and vacuum polarization graphs are built up in QED. 
The counterterm is in the image of R, while </>(r) -|- Yl-^izv ^'Ri'^)4'(^ H) is the result of 



Bogoliubov's famous R(T) operation on T which eliminates subdivergences in F [13|. 

Under suitable conditions on the behavior of ii, in this ratio of characters short distance 
singularities drop out [U, |l^ . 

Before we comment on the renormalization scheme chosen for our calculations, let us 
introduce the bidegree of a Feynman graph. This standard notion can be introduced for any 
Hopf algebra which is reduced to scalars ||^, |^ by the counit e with 

^(^l) = q-, e(^) = 0, else. 

If we decompose H = Hq (B -ffaug, with -ffaug being the augmentation ideal as the kernel of e 
we can investigate, for any positive integer k and Hopf algebra element X, 



Xk := A''-\X)nH, 



aug- 
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For sufficiently large k this will necessarily vanish. We define the bidegree bid(X) as the 
largest k such that Xk 7^ 0. Elements in Hq have bidegree zero. Note that Hopf algebra 
elements of unit bidegree are precisely the primitive elements X € -ffaug in the Hopf algebra, 
with 

A(X) =X®l + l(g)X ^ i^aug ® ^^aug- 

Having introduced this standard notion we introduce the renormalization scheme for which 
we choose minimal subtraction (MS). Each application of the scheme is given by a projection 
onto the pole-part of the considered Laurent series and is symbolized with brackets "()": 

+00 -1 

j:=-r j-=~r 

The reader should convince himself that this map makes the ring of Laurent series with poles 
of finite order into a Baxter algebra, 

{ab) + {a){b) = {a{b)) + {{a)b). 

Note that the degree r of the pole terms in a Laurent series assigned to a graph by the Feynman 
rules in dimensional regularization is in general majorized by the bidegree r < bid(r) and 
equals the bidegree in our simple applications. 

We expect to encounter Cl'^) inside the coefficients of the Laurent series in the regulariza- 
tion parameter e emerging from a series expansion of the functions -Fa,6- 

In QED we have to take our matrix-calculus into account. The renormalization of such 
matrix expressions is now given by inserting a diagonal matrix R, which consists of the 
renormalization map as entries: 

V U Kms J 

Inside a string of matrices, this matrix has to be inserted wherever the renormalization map 
is applied. It acts on expressions on the right. 



4 Rooted Trees 

Before we build up graphs and calculate their counterterms, ie. their antipodes, let us first 
mention that in our simplified context, the Hopf algebra of Feynman graphs is isomorphic 
to a Hopf algebra of rooted trees with a very small set of decorations given by our one-loop 
graphs. 

In the class of graphs to which we have restricted ourselves one-particle irreducible sub- 
graphs are either nested in each other, or disjoint. They hence form tree-like hierarchies. 



and one easily translates graphs in rooted trees |15], with a one-to-one correspondence be- 
tween one-particle irreducible subgraphs and vertices in the rooted tree (the map from graphs 
to rooted trees is one-to-many for overlapping divergent graphs, and can be systematically 
constructed [jl6|): 

The translation from a Feynman diagram to a rooted tree has to be done in the following 
way: Set a box around the subdivergences of a Feynman graph T and mark the upper hor- 
izontal hne with a dot (~ vertex). Dots of nested boxes, that is boxes where one of them is 
contained inside the other, are connected with a hne (~ edge) (see fig. P). 
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Figure 5: Translation of Feynman diagrams into rooted trees. Here we are dealing with decorated 
rooted trees, where the subdivergences are written next to the vertex to which they belong. In the last 
column one can see another way to encode the structure of the trees in the form of nested lists [ p4| : 
The entries are again the divergent parts of a Feynman diagram and the formation of the parentheses 
gives the structure of the tree, beginning with the vertices "at the end" of a tree up to the root. 



Each vertex of such a tree represents a divergent subgraph of the diagram. The antipode 
of such a tree is a sum over all full cuts or, recursively, iterates over admissible cuts using the 
coproduct. This easy-to-implement Q coproduct 

A(r) =r® i + i0r+ ^ p"(t) i2^(T) 

admissible cuts c 

is by now a standard combinatorial tool (see, for example ||l^). Using it, the counterterm is 
a recursive construct 

Sms{T) = -(r(r) + Yl SMs{P'iT))R'iT)), 

admissible cuts c 

where an evaluation of the graph corresponding to the tree using the Feynman rules is un- 
derstood before applying the projection () for minimal subtraction (MS), in accordance with 
the general formula given for graphs. 



5 Implementation 

We shortly describe an implementation of our scheme in a computer program. This imple- 
mentation is not entirely self-contained in the sense that given a Feynman diagram, it would 
compute the antipode. Instead, its input is a decorated rooted tree in list notation such as 
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(r(n)(E)). The actual construction of that rooted tree from the diagram, as described in 
section |^ is left to the user (see also fig. P). 

Our implementation uses the GiNaC system for symbolic computation in the C++ pro- 
gramming language . GiNaC provides efficient implementations for handling Laurent series 
as needed in dimensional regularization. From C++, all the linguistic instruments for object- 
oriented programming are borrowed and available to us. We follow a traditional approach for 
representation of our rooted trees where there exists a container class called node that may 
or may not have several children of the same class. In addition, each vertex in an object of 
class node contains an object of the abstract and polymorphic class decoration. From the 
abstract decoration class a number of concrete classes like Sigma, Gamma and Vacuum are 
derived, corresponding to primitive Feynman diagrams S, F and 11, as well as some additional 
classes that allow us to distinguish between the position of a subgraph inside its supergraph. 

The trees under consideration are not ordered with respect to their children, i.e. the trees 
(r(E)(n)) and (F(n)(S)) are equivalent. Therefore the children form a multiset where only 
the multiplicity of occurrences is relevant. 

A template class multiset is part of the C++ Standard Template Library (STL). It has 
the additional advantage that the elements are always automatically sorted with respect to 
some specified ordering. This turns out to be useful for convenient identification of equivalent 
nodes and also to establish an order relation on them. The state of the edges leading to the 
children of a node (either "cut" or "uncut") needs to be taken into account as well, so the 
multiset is really one of pairs of nodes and boolean variables. We chose to delegate methods 
from class node to the corresponding class derived from decoration using dynamic type- 
dispatch. Hence, the decoration must be stored as a pointer, calling for some hand-made 
memory management in class node. A completely realistic layout of our class node is then: 

1 class node { 

2 public : 

3 // constructors, destructors, delegators, 

4 // etc . . . 

5 private : 

6 decoration *deco; 

7 niultiset< pair< node, bool > > children; 

B }; 

Note that a node can be either the entire tree, or a subtree or a single (atomic) leaf. The 
layout of class decoration holds a pair of indices and in the case of QED a GiNaC expression 
(class ex) for the gauge. An object of class ex is entirely sufficient, since it may contain 
either a numerical value (integer or fractional) or a symbol (like .^) (or even more complex 
expressions, if need should arise). 

The knowledge how to manipulate the indices i and j, depending on the type of deco- 
ration and on the state of the edges is built into the classes of type decoration. They are 
automatically adjusted when the tree is constructed. Furthermore, the state of the edges is 
also irrelevant for the user of the program since we need all possible combinations of cuts. 
If there are n vertices inside a tree, then there are n — 1 edges and we have to construct all 
possible 2"~^ combinations. 

Once all the trees have been created we call a method called evaluate in each of them that 
traverses the tree in a top-down fashion: In each node, the expression in terms of Fa-\-ie,b+j£ 
is dispatched and expanded as a series in the regularization parameter e before evaluate is 
called on the multiset of children. The resulting 2"~^ Laurent series are then added together 
and the coefficients are expanded. This procedure yields the antipode. 
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Here is an example how the programs are used in practice: 



1 $ ./qedl " (((Sigma [0] ) (Sigma [xi] ) Gamma [-1] ) Gamma [xi] ) " 

2 After decoration the tree has these indices: 

3 (Gamma [3 , 0] [xi] (Gamma [2,0] [-1] (Sigma [0 , 0] [0] ) (Sigma [0 , 0] [xi] ) ) ) 

4 + # + # + # + # + # + # 

5 The antipode of this tree is: 

6 

7 (-l/4-l/4*xi-2-l/2*xi)*x-(-3) 

8 +(7/24+l/24*xi-2+l/3*xi) *x- (-2) 

9 +(-5/16-13/48*xi-2-7/12*xi)*x-(-l) 

The graph to be computed is passed in tree form as a string on the command line, together 
with the gauge parameter. Note how the indices i and j are then set up automatically in 
line 3. The next line is a simple progress bar, useful when computations become longer. The 
result is then printed as a power series in the regularization parameter, called x instead of e 
in order to please the computer. 

Computationally, the results are always quite small (c.f. section but there is a huge 
swell of intermediate expressions. After all, the Laurent series arise from expanding products 
and sums of lots of F-functions and so inevitably carry Euler-Mascheroni constants. But these 
have to vanish by cancellations, a fact that is conveniently used as a consistency check for the 
result. In practice, antipodes of nine-loop graphs in massless QED are still computationally 
feasible but may require several hours and hundreds of megabytes memory. This emphasizes 
the drastic loss of efficiency in comparison with |^, ^, which is compensated by a gain in 
flexibility to handle different decorations at any step in the recursion. 

The programs developed herein are written in plain ISO C-|--|- and will run on any 
system where the GiNaC library has been ported to. Note that porting to a new platform 



will require porting CLN |18] to that platform first, since GiNaC depends on CLN for its 



arbitrary precision arithmetic. 



6 Computational Results 

In a first step, to examine the appearance of ^-functions, we will now calculate the antipode 
of different graphs given in figure ^ and figure |^ that represent the two topologies of figure |l|. 
In Yukawa theory the graphs T^^ ^ and T'^^ ^ ^ of figure ^ are representations of the ladder 
topology l^.b, the graphs q and ^ -^j of the swiss cheese topology |l].a. The same holds 
for QED with [1, 1, 1] replaced by [1, 1, 2] and [2, 0, 1] by [2, 0, 2]. 

The graphs for the Yukawa theory, together with their expressions in characterizing func- 
tions are given in figure ^. The figure also contains the decorated rooted trees we get from 
them to calculate the antipode. One can see that except from the decorations they all belong 
to the same rooted tree, and hence have the same structure as far as their short distance 
singularities go, as advertised in the introduction. 

In QED, the considered graphs in figure ^ possess one more vertex correction compared to 
the graphs in yukawa theory. This additional photon line is necessary because the short dis- 
tance singularity structure is actually determined by the bidegree, ie. the number of subdiver- 
gent graphs. It so happens in QED that the one-loop fermion self-energy and one-loop vertex 
corrections (with divergent subgraphs) are overall finite if the internal photon is transversal. 
In our case, we thus find that the insertion of a one-loop vacuum polarization into the internal 
photon line in a vertex results in a convergent vertex correction. We have to plug the whole 
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Figure 6: Graphs in Yukawa theory and their decorated rooted trees. 



function into one more vertex correction to get to the next level in the bidegree. And indeed, 
we need an additional vertex correction, compared to Yukawa theory, to obtain the ("(3) in 



the pole-terms: gauge symmetry delays the appearance of transcendent als |]19[. Also let us 
mention that we can easily compare graphs which have an internal vacuum polarization with 
graphs which have an internal fermion self-energy by using the before-mentioned elimination 
of the trace in the vacuum polarization, using 



no,o = tr(l)iFi,i 



tr(l)So,o. 



Furthermore, F^-j^ ^ 2] 2] ™ QED do not have this extra shift between loop number and 

bidegree by themselves, as there the desire to maintain the same topology never forces us to 
use a vacuum polarization as a subgraph: all one-loop subdivergences are fermion self-energies 
and vertex corrections. However, we use the same gauges in these graphs as in the graphs 
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Figure 7: Graphs in QED and their decorated rooted trees. 

Fj^j^ ^ 2] ^-iid Fj2 Q 2] to compare the results directly (see below), and therefore obtain the same 
difference between loop-number and bidegree. 

Finally, in our results, we are only interested in residues of counterterms, poles of first 
order. All interesting relations between transcendental degree and topology will appear there. 
The scattering type formula Q will make sure that parts of these relations will then resurface 
in the poles of higher order in the counterterm, but they contain no new information. So 
in the following we solely exhibit the residues of the counterterms for our selected class of 
graphs. We denote by res(r) this coefficient of the pole of first order in the MS counterterm 
in dimensional regularization of a graph F. Assorted results for these residues are: 
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Figure 8: Iterated vertex corrections. 



Yukawa theory: 

.2 ^ 
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QED: 

^4,1.2]) = -^4^(301 + 143^-170^2)^ 

Mrf,,i,2]) = -(^(584 + 659^ + 59^2)) 

res(ri2 o_2]) = - (^(521 + 309^- 236^2) ^^^(3)(^^^)2) 

res(r|o_2j) = -^-L(832 + 879e + 3ie2) + ^C(3)(l + e)') 

We see, that in Yukawa theory the graphs q ^ and q ^ with the swiss cheese topology 
have a residue involving C(3), while Fj^^ -j^ and V'^-^ ^ realizing the ladder topology just have 
a rational residue, as expected. A similar result holds in QED: the graphs q 2] ^^'^ ^ 2] 
with the Swiss cheese topology again have a residue involving C(3), while T^^ ^ 2] ^fi 1 2] 
realizing the ladder topology just have a rational residue. Those residues are in general a linear 
combination of terms of varying transcendental weight. This weight vanishes for a rational 
number, and in accordance with standard practice in the study of multiple zeta values the 
transcendental weight of a monomial CiJi) is 

i i 

Then, the above results confirm that the coefficient of the highest weight transcendental in 
the transition from res(rj2 q 1]) to res(rp ^ (res(rj2 q 2]) to res(rp ^ 2]) for QED) is invariant. 
This is the partial symmetry announced in the introduction. It holds for any gauge in the 
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QED case and holds in general, if we increase the loop number as for example in figure (^). 
In Yukawa theory we get: 

= - las) - |;C(4)) 

-ff|,„„) = -(l + i^cp) + i^a4) + ^c(5) + 4c(3)^ + ^«6)) 

/^i X /^6481 33613,, , 2133,,,, 101 , , 27 ,n 

27 3 7 \ 

+— C(6) + ^C(3)C(4) + ^C(7) 
896^ ' 4480^ ' 1920^ 7 

Similarly, for QED we find: 

res(rf2o_2]) = -^(521 + 309e- 236e')-^(l+O'C(3) 
res(rf2^0 2]) = -^(832 + 879e + 31^2) -^(l + e)'C(3) 

-|(i + em4) 

res(rf2,o,3]) = -^(^3327 + 32578^ + 26038^' + 6811^') - ^(1 + 0'(31 + 40C(3) 

-|^(i + em4) 
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res(rf2 4]) = ~4o^(^^''^^^ + 1048812^ + 889375^^ + 130774^^ - 145630^^) 

~8i0^^ + ^^'^^^^^ + ^^^^ + 736C')C(3) + 2^(1 + 0'(-121 + 230C(4) 
-^(1 + 0^C(5) 

^^^(^[2,0,4]) = - ^gj^^gQ (1357764 + 4016396^ + 4028301^^ + 1438370^^ + 68749^^) 

^ -(1 + 0^(1448 + 23654 + 836S,'^)Ci3) - -^(1 + 6^(29 + 20C(4) 



1680' ' " - /^v / ;L40 

-4(l+«'«5) 

^es(rf2,o,5]) = 430295040 (-19150607852 - 55087254529^ - 55376777218^^ - 24306202368^^ 

-10014576786^'' - 5440770359^^) 
+^(1 + e)'(-39057 - 40975^ - 19389^2 ^ 6061^3)^(3^ 

^ :i + 0^(2031 + lie + 446^2)^(4) + + 4)4(_347 + 670C(5) 



4480' ' " " ' 960 

^ (i+e)W- 7^7(1+ erne) 



1120' ■ ^ 224 

1 

1290240 



^^^(^[2,0,5]) = - ,»Q^^.^ (29658556 + 111493999^ + 1615393734^ + 1128129934^ 



+386883794^ + 55793644^) 



(1 + 4)^(80770 + 1555504 + 795214^ + 34724^)C(3) 



26880 

-^(1 + 4)'(3254 + 39194 + 13404')C(4) 

-^(1+ 4)^(179 + 84)C(5) 

^ -(1 + 4)5^(3)2 -A(i + ^)5^(6) 



1120' ^' " ' 224 

The reader will note that in QED our residues are polynomials in the gauge parameter of 
a degree reduced by two steps from what one might expect, to enable comparison between 
configurations with insertions of self-energies into photon or fcrmion lines. The corresponding 
fermion self-energy was for that purpose evaluated in the Fcynman gauge, and in the Lan- 
dau gauge for the affected photon propagator. If one only compares cases with self-energy 
insertions at fermionic lines, one can abandon these restrictions and we did confirm that the 
reported invariance holds as expected with coefficients which are polynomials in the gauge 
parameter of degree equal to the number of photon lines. 
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7 Discussion and Proof 



Actually, results of the form reported in the previous sections can be derived from the analytic 
structure of the functions Fa^i, and some basic field theoretic arguments. 

Let us reconsider the situation. The simplicity of Feynman graphs considered here mani- 
fests itself computationally by the fact that they factorize in a unique manner. Each divergent 
subgraph 7 depends only on a single external momentum q say (that is the reason why we 
only consider vertex subgraphs at zmt), such that its evaluation in dimensional regularization 
gives a result of the form 

i 

Here, ^(7) is the number of loops in 7, and Fi{e) are g-independent form factors, and the 
dimensionless Ci{q) are 

= C2(g) = q^i/q^ for the QED vertex (the only case here in which the sum has more 
than one term), 

ci{q) = 1 for the vertex correction in Yukawa theory, 
ci{q) = <ji for any fermion self-energy, 
ci(9) = - ItiQu/q^ for the photon, 
ci(g) = 1 for the scalar boson. 

Insertions of such graphs 7 in another graph T only raises powers of the scalar part of 
some propagator of F: 

1 1 

^ ^ [g2]l+n(7)e- 

We can keep track of this by notating these loop numbers n{"f) in the entries M(F).y for the 
corresponding propagator in the adjacency matrix M(F). 

Let us now assume that F is some primitive vertex correction, ie. free of divergent sub- 
graphs, and let us write as before F^ and F^ for two distinct choices of zmt. 

Consider a bunch of IPI graphs X = Y^=i 7i each dependent on a sole external momentum 
as described above. Let k = h\d{X) be the bidegree of X, so that X has a highest pole in e 
of degree k with coefficient c^. Let now Gx be chosen gluing data such that Fx := ^ '*^Gx 
is obtained from inserting X at specified vertices and propagators into F, with bid(F) = 1 
without loss of generality. (Any IPI graph can be written in the form T-kQ^ X for appropriate 
such F, X ||l^, |l^, in generalization of the closed Hochschild one cocycle i?+ of undecorated 
rooted trees). Further, each X allows for an expansion 

ct>{X) = i,{l + T{X){e)) 

and similarly, let (/'(F) = £2^0.(1 -|- T(F)(e)). Now assume that the Taylor series 
[1 T{X){e)\\V T(F)(e)] = 1 EjliCj-e^ is such that the transcendental weight w{cj) 
of Cj increases with j: 

w{cj) < w{Cj+l), 

Vj > k — 1. Here, we define the transcendental weight of an expression which is a sum of 
terms as the highest transcendental weight appearing in its terms.0 

''The question as to how define the transcendental weight in a context which exceeds the Riemann C^- 
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Proposition 1 The counterterm is the same for and F^, and hence their residues are 
equal. 

Proof: Elementary, as F^ — F^ is UV convergent, and hence F^ and F^ generate the same 
counterterm. □ 

In particular, we also note that in the above, F^ — F^, when inserted into another graph, 
produces a result with a bidegree reduced by one unit compared to the insertion of either F^ 
or F^ alone. 

This is not yet the desired result, as in our case we have to compare F^^ with F^^, where 
Xi is a collection of subgraphs in which all vertex subgraphs are of type F^, and X2 is the 
same collection of subgraphs apart from the replacement F^ — > F^ for all vertex subgraphs. 

Any graph of type F^^ or F^^ , which itself can contain subgraphs Xi of these varying types 
of vertex corrections plus self-energy subgraphs, can now be expressed in terms of the other. 
Similarly, this holds for these vertex subgraphs of either type, on the expense of generating 
extra terms of reduced bidegree 

bid(Fi,^ - r2,J < bid(Fi,J = bid(F2,J, 

which involve differences F^. — F^. for appropriate Xi. Hence, under the above assumption 
of monotonic increase of the transcendental weight with the bidegree, we get upon iterating 
such insertions 

Proposition 2 The coefficient of the term of maximal transcendental weight is the same for 
res(F^^) and res(F^^). 

Here, Xi,X2 are related, as above. 

This explains immediately our results as a look at the functions Fa^b, and hence the 
corresponding evaluations of our subdivergent graphs, shows that they fulfill the required 
assumptions of monotonic increase of transcendental weight, which was completely determined 
from the appearances of the Riemann (^-function in our simple examples. Note that the 
factorizations into two-point functions and the absence of all other primitive graphs apart 
from one-loop functions were the two main simplifications which enabled us to satisfy the 
assumption. 

The study to what extent a sensible transcendental weight can be established in general 
will be a topic of future work. Any sensible answer we will expect to deliver the same 
permutation invariance of the residue as reported here. 
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functions or MZVs we do not have to answer here. Also, the attentive reader might have noticed that 
we set the transcendental weight of the gauge parameter to be zero for the QED results, treating it as an 
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